CLOSED FORMS FOR A MULTIGRAPH ENUMERATION 



CHRISTIAN BROUDER, WILLIAM J. KEITH, AND ANGELA MESTRE 



Abstract. We partially sum a recursion formula related to the Hopf alge- 
braic generation of connected Feynman graphs in quantum field theory. One 
of these sums is Sloane's sequence A001865 in The On-Line Encyclopedia of 
Integer Sequences. This formula is a particular case of a more general recur- 
sion formula for a function counting inverses of automorphism group orders 
for multigraphs. 



1. Introduction 

Denote by I{n, k) the sum of the inverses of the orders of the groups of auto- 
morphisms of ah the pairwise non-isomorphic connected muhigraphs on n vertices 
and cyclomatic number k. The third author, in joint work with Oeckl [3], used the 
symmetric algebra on a vector space to derive a recursion formula for generating 
certain linear combinations of tensors over the rational numbers. The tensors rep- 
resent connected multigraphs (with loops and multiple edges allowed) on the same 
number of vertices and on the same cyclomatic number. The formula produces 
larger multigraphs from smaller ones by increasing by 1 the number of their ver- 
tices or the number of their edges. The key feature is that the sum of the coefficients 
of all the tensors representing isomorphic multigraphs is the inverse of the order of 
their group of automorphisms. Note that, in this setting, the order of the group 
of automorphisms of a multigraph consisting of an isolated vertex with s loops is 
2*s!. Also, the order of the group of automorphisms of a multigraph consisting of 
two vertices joined by r multiple edges is 2r\. 

An interesting property of the formula of [3 is that of having an alternative 
expression: see Proposition 15 of that paper. The latter consists of two terms. The 
first one relates connected multigraphs on n vertices and cyclomatic number k with 
all their connected subgraphs on n vertices and cyclomatic number fc — 1. The 
second term relates connected multigraphs on n vertices and cyclomatic number k 
with all the pairs of their connected subgraphs with total number of vertices equal 
to n and total cyclomatic number equal to k. This alternative formula then induces 
a recurrence for /(n, k). 

This note proves closed or at least simpler forms for the sequences of I{n, k) in 
constant n or fc when these are small. In constant fc, these are 

Theorem 1. For all n>l, I{n, 0) = and I{n, 1) = 5 ELi 
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The first formula is the sum of the inverses of the orders of the groups of auto- 
morphisms of all the trees (up to isomorphism) on n vertices, which was pointed 
out by Polya [5]. 

In constant n, the answer is much more complete. We find several forms of the 
generating function 

oo 

Rnix) = ^/(r^,/c)a;^ 

including the two- variable generating function 

Theorem 2. ^^=1 = a; log (^^=0 7^e™V2) _ 

Define J{n, k) := 2''{n+k — l)ll{n, k) to normalize the sequences I{n, k). We find 
that the sequences of J in constant n possess concise rational generating functions 
whose coefficients are simple formulas. With the sum running over all compositions 
of n, i.e. sequences of positive integers that sum to n, the general form is: 



Theorem 3. 



1 " 



i-iy 



There are 2"~^ compositions of n, but the number of distinct sums of squares 
of parts among these compositions is much smaller, meaning that the number of 
factors in each generating function is much smaller than might have been expected 
from the above theorem. Hence a more efficient method of producing a generating 
function and explicit formula for a given n might be desired. There is a fairly simple 
routine for calculating explicit formulas and generating functions, which we give in 
the last section. Values for n = 1 through 4 are given in the table in the theorem 
below: 

Theorem 4. For k > 0, the formulas and generating functions for J(n,k), 1 < 
n < A, are as given in the following table. 



J{n, k) 


Formula 


Generating function 


J(l,fc) 


1 


Sfc=o '^(1' ^)^'' ^ 


J(2,fc) 




1/2 

(l-22;)(l-4a;) 


J(3,fc) 




(l-32;)(l-52;)(l-92;) 


J(4,fc) 


(f ) 16'= + (f ) 6*^ - (^) 10'= - (2)4'= - (8)8'= 




(l-4a;){l-6a;)(l-8£c){l-10a;)(l-16a:) 



In the next section of this note, we prove Theorem 1. When fc = 0, we use 
Dziobek's recurrence for Cayley's formula to show that J(n,0) — n"~^. When 
fc = 1, we use an identity related to Abel's theorem to prove that J(n, 1) = 
"-'Sjl^i "n-^)! 1 which is Sloane's sequence A001865 in |5j. In Section |3] we prove 
Theorem 2 by solving the differential recursive equation. We also give a general ex- 
pression for J(n, fc). In Section |4] we consider i?„ from another angle, and re-prove 
some of the expansions with a more hands-on approach that displays different un- 
derlying combinatorics. 
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2. Formulas in constant k 

We use the formula of Proposition 15 of ^ to give a recursive definition for our 
sequences, which in Theorem 1 are respectively A007830 and A001865 in [S]. 

In the following, for all n > 1 and A: > by I{n, k) we denote the sum of 
the inverses of the orders of the groups of automorphisms of all the pairwise non- 
isomorphic connected multigraphs on n vertices and cyclomatic number k. By 
the formula of Proposition 15 of [3] the numbers /(n, k) are defined recursively as 
follows. With boundary condition 

/(n, fc) = for fc < and/or n < 1 , and /(I, 0) = 1 

then for all n > 1 and fc > we have 

(1) /(n,fc)= ^ 



2{n + k- 1) 

n-1 k 

I{n, k -1) + ^^i{n- i)I{i,j)I{n ~i,k-j) 
1=1 j=o 

For convenience of calculation, we define 

(2) J(n,fc) :=2'=(n + fc-l)!/(n,fc). 
Then, for all n > 1 and fc > the recursion for J{n, k) is 

(3) J(n,fc) = nV(n,fc - 1)+ 

n — 1 k 



1=1 j=0 

The same boundary conditions hold. 

We now prove Theorem 1. 
Theorem 1 (restated). For all n > 1: 
(a) J(n,0) = n"~^ , or equivalently, l{n,0) = — 



1 fn + k ~2\ 

2 ( - I ■_i ]^in-i)Jihj)Jin-i,k- j) ^ 



(b) J (n, 1) = nl J2l=i ' or equivalently, lin, 1)^2 TZ=i JJI^. ■ 

Proof. To prove formula (a) note that the recurrence for J(n, 0) reads as follows: 

J(1,0) = 1, 

-^^""'^^ ^ 2n{n - 1) E ( J^'(" - »)V(z,0) J(n - », 0) . 



i=l 



For T{n) :— nJ{n, 0) we obtain Dziobek's recurrence [T| for Cayley's formula: 

r(i) = 1, 

1 ""^ fn\ 

(5) T{n) = (^Jj(„-z)T(*)r(n-*). 



4 



BROUDER, KEITH, AND MESTRE 



This is thus solved by T{n) = n"-^ which is the sequence A000272 in [5]. There- 
fore, we have J{ti, 0) — n"^^. 

To prove (b) we first recall the following Abel- type identity [4j p. 93]: 

n{n + y)"-^ = ^ ("^^ i{-x + l){~x + iy-^{x + y + n- 
i=i 

For a; = and y — j the above formula specializes to 



n 
i 

We now proceed by induction on n. The result clearly holds for n — 1. We 
assume the result to hold for alH in < t < n. Then 

^<M, ^ .-..i|:((»:;).-,„-.„|:tt. 



n - -^z! > - 



2n -^-^ V i 

2—1 



n ^ \ i J (n — i — uV. 



n — ' — ' t\ in — I — u)\ 

^—1 i—l 

Now, formula ^ yields for n = m + ^: 

^-^ i'. (m — i)\ ml ^-^ \ i J ml 

Substituting this into the previous line we obtain 



J(n,l) = n"-^ + in!^ ("7^)"7' =.»-^+n!^ 



(n-^)! ("~/^-l)■ 



V — = n' V 



as desired. 



□ 
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A natural question is whether the recurrence can be similarly simplified for fc > 1. 
Attempts with Abel-type identities of greater generality did not yield results that 
we considered useful, but the quest does not seem theoretically implausible; the 
authors would be interested in any observations that readers might be able to 
provide. 



3. Formulas in constant 7i 

3.1. Generating functions. Define the generating function for the sequence of 
I{n, k) with n constant: 

oo 

(7) Rn{x) = ^/(n,fc)x^ 
In particular, 

(8) = f; /(I, k)x'^ = £ ^^x'^ = e-/^ 

k=0 k=0 

because J(l, k) — 1 for all k. 
We now wish to prove 

Theorem 2 (restated). 

Rn{x){sxy' = xlog ^ 



n— 1 \n— 

Proof. If we multiply each /(n, k) by x^ and sum, the recursive equation ([T]) for 
/(n, k) becomes a differential recursive equation for i?„: 

ji-i 

2{n - l)Rn{x) + 2xR'^{x) = n^xR^ix) + Yi{n ~ i)Ri{x)Rn-t{x), 

i=l 

with boundary conditions Rn{0) — I{n, 0) — n"~'^/nl, Rn{x) = for n < 1. 
To simplify, apply the change of variable 

pnx/2 

The inverse relation 

= n!(logz)"-iz-"/2^„(logz) 

gives the boundary condition Pri(l) = Sn,i- 

The differential recursive equation satisfied by Pn{z) is 

n-1 y s 

2zP;(z) ^ n{n-l)Pn{z) + Y,Kn-^)^^P^{z)Pn-^{z). 
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Observe that Pi{z) — 1, and if Pi{z) is polynomial for all 1 < z < n, then P„(z) 
satisfying this recurrence is polynomial as well. The first polynomials are 



Pi{z 

P2{Z 

Psiz 

P4{Z 
P5{Z 

Peiz 

P7{Z 



= 1, 

= -l + z, 

= -6+12z-5z^ -4z^ + z*^, 

= 24 - 600 + 30^^ + 20z^ - lOz^ - 5z^ + z^°, 

= -120 + 3600 - 2700^ - 900^ + 120^^ + 20^^ - 15^^ - 6z^° + z^^, 

= 720 - 25200 + 25200^ + 2100^ - 12600^ + 2100^ - 700^ + 2100^ - 350^ 

+420^0 -2I0"- 70^5 + 0^1. 

The equation for the generating function 



/(y,0) = yPniz) 



n=l 



r 



IS 



2z^ - 2^ ( 

dz ^ dy"^ \ dy) ' 

with the boundary condition 



= Y^P^ilf =y. 

ri ~ 1 

The non- linear term is eliminated by the change of variable / = loggi. The 
equation for g is therefore linear: 

dz dy^ ' 

with the boundary condition g(%), 1) = e-^^^'^^ = e^. A family of solutions of this 
equation is given by 

n\ • 

Moreover, the boundary condition g(y, 1) = implies gn = ^ for all n: 

yn^n{n—l)/2 



Finally, 



n=0 



f{y,z) = E^"W|!=l°g(E 

n—l 

Note that if we set = (1 + f), we obtain 



°° yn^n{n^l)/2 



n=0 



/ °° „«n _|_ f^"(n-l)/2 s °° ,,n 

/fe. ') = >»s ( E ) - E onw^. 

71=0 ■ 71=0 
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which is the well-known exponential generating function for Cn{t) = '^Qt'^^'^\ 
where the sum is over all connected graphs G on the set {1,2,..., n} and e{G) is 
the number of edges of G. (See [51 for this formula and related results.) Thus, 

Pn{z) = Cn{z-l). 

Substituting for z, we obtain a generating function for R^'. 



OO ^-j^ ^ CO 

^-^ nl X ^-^ 

n— 1 n— 1 

The generating function for Rn{x) becomes 



-x/2\n 



J2Rn{x)s- = :,/(se^/Vx,e-) = xlog(^if^e™^/2 

n—1 n—0 

However, this expression is highly singular in x and does not provide a generating 
function for Rn- The alternative expression 



OO OO 

ch{x,s) = ^i?„(a;)(sa;)" = 2;/(se^/2^e^) a:log(^ 



n=0 



/2 

n' 



is not singular and generates in the sense that Rn{x) is x " times the coefficient 
of s" in the expansion of xh{x,s). □ 

The extraction described yields 



m — ' ni! . . . ri„ 

m=l rii4 |-nm=n 



where all n, > 1. 



3.2. Explicit forms for J{n,k). We now wish to show Theorem 3 by isolating 
explicit generating functions for J{n, k) individually. We can proceed as follows: 



OO OO 



xh{x,s) = /(n, k)s 

n=l k=0 

Now replace x — >■ 2x, s — >■ s/2 to obtain 



OO OO 

+fc 



2xh{2x,s/2) = ^ ^ 2'=/(n, fc)s"a;" 

n=l fe=0 

We define /3(s, t) = 2 dxe-''h{2xt,s/2) to get 



nOO ^ ^ /"OO 

2/ dxe-''h{-2xt,s/2) = J2J2 dxe~'=2''I{n,k)s"{-xt)"+'"^ 



/3(s,-t) 



J2J2'^''{n + k-iy.iin,k)s^\-ty 



\n~\-k—l 

where we used e~^x^dx — p\. Therefore, 



OO OO 



n=l k=0 
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is a generating function for J(n, k) and this equation defines Zn{t) — X^fe^o "^("'' 
We calculate 



/3(s,-i) 



•^0 n=0 ' n=0 



where we have put A = e ^ . This gives 

" (-1)™ ^ (./2) 



"iH bn„ 



m ^ — ' ni! . . . n„ 

m— 1 ni,...,nT7i 



Thus, 



m ^ ni!...n„! 1 - (n^ H ^ n^Y 



The polynomials Zn(t) can be calculated by summing over the compositions of 
n, proving Theorem 3: 

(-1)™ 1 1 



(2t)"-i ^ m ni!...n„!l-(n2 + .-. + n^)r 

^ ^ m— 1 niH I-tIth— n ^ ^ 

□ 

The number of distinct sums of squares of the parts of partitions or compositions 
of n is the OEIS sequence A069999; it grows like n? /2 [7]. (We remark that the 
first case where a sum appears twice is n = 6.) The first Zn{t) = X]fe>o "^("•' ^)^'^ 
are: 



Zi{t) 
Z2{t) 

Zi{t) 

z^it) 

Zeit) 



1 



1-t' 

1/2 



(l-2<)(l-4t)' 
1 

(1 - 3t){l - 5t){l - 9t)' 
4-34i 



(1 - 4<)(1 - 6t){l - 8<)(1 - 10i)(l - 16t) ' 

25 - 606t + 3557^2 
(1 - 5t)(l - 7t){l - 9t)(l - - 13t)(l - 17t)(l - 25i) ' 

24(9 - 451< + 7292t2 - 37m0t^) 
(1 - 6t){l - 8i)(l - 10<)(1 - 12i)(l - 14i)(l - 18t)(l - 20t)(l - 26i)(l - 36i) ' 



4. Formulas in constant n: the hands-on approach 

In this section we are interested in re-deriving the Zn(t) more directly, and 
finding the closed formulas to prove Theorem 4. Using strictly discrete mathematics 
instead of the generating function, we seek additional insight into the combinatorial 
structure of the sequences J(n, k). 

The case J(l, fc) = 1 for fc > is trivial from the recursion. 

For n — 2, the recursion (jSj becomes 
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J(2,fc) = 4J(2,fc-l) + i^ ('':]j(l,j)J(l,fc-j-) 
= 4J(2, k-l) + ^Y.(^) = 4^(2' + 

3=0 

We now note that 4(4'=-i - 2'=-2) + 2'="! = 4*^ - 2'= + 2'="! = 4*= - 2'''^i and since 
the statement holds for fc = 0, the claim on the values follows inductively. 

That the value 4^^ — 2''~-^ is the coefRcient of x'^ in 7771 — tt-th — can be seen from 

2(1 — 2a:)(l— 4x} 

expansion: the coefficient is 



1 ^^k _,_ 4fe-i2 _j„ 4/C-222 








1 

2^ 


1 

4 + - 


1 

■ + ^ 










2^, 


1 =4'= 





For n = 3, the recursion (jS]) gives 



(9) J(3,/c) = 9J(3,fc - 1)+ 



fc 

= 9J(3,fc-l) + ^ 



^ + ^) J(l, j) J(2, - + Q ^ J ) J(2, j) J(l, fc - j) 



3=0 



k + 1 



fc + 1 



(4^ - 2^- 

=.9J(3,fc-l) + 25:(';:M(4^-2-i) 

1=0 ^ 



2 — \ I 

i=0 



fc + 1 



9J(3, fc - 1) + i(l + 4)'=+i - ^(1 + 2)'=+^ 



9J(3,fc-l) + i5^+i-i3'=+i. 



Let us pause to note the process here. We added a term to the sum, namely the 
i = term, which is the j = —1 case of 4^ — 2^^^. However, this term is 0. We 
then complete the binomial sum, changing the 4 to a 5, and the 2 to a 3. 

Iterating the recursion with the form thus obtained, we obtain the closed formula 
for J(3,fc): 
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J{3,k) 



9J(3,fc- 1) 



-k+1 



ifc+1 



9V(3, - 2) + 9 • i (5^' - 3^) + i (S'^+I - 3'^+^) 

• • ■ = 9^= J(3, 0) + i (g'^-i (52 - 32) + ... + (5^+1 - 3^=+!)) 

igfc+i . Q (5I _ 3I) + . . . + 9-(fc+i) (5fe+i _ 3^+1) j 



2 



Tc+l . 



To show that the generating function claim is true, we factor the last term of 
equation ([9]): 



J(3,fc) 



9J(3,A:-l) + i5'=+i-i3'=+i 
1 

9J(3, - 1) + ^(5 - 3) ^ 5''-'3' 



i=0 



9J(3,fc- 1) +^5'=-*3* 



i=0 



The same recurrence and initial condition is satisfied by the coefficient of x'^ 
in (i_3a;)(i_\;j)(i_9a;) J we can see by noting that this coefficient is the complete 
homogeneous symmetric polynomial in 3, 5, and 9: 



L^J(l~3x)(f-5x)(f-9x) 



9 J2 9'^5''3" + ^ 5'=-*3' 



i=0 



where [x'^] f{x) denotes the coefficient of x*^ in f{x). 

The n = 4 case is proved similarly. Begin by expanding the defining recursion for 
J(n, k) using the formulas for J(l, k) through J(3, k). Gather symmetric expansions 
with the same binomial coefficients. 
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J(4,fc) = 16J(4,fc- 1) + 




Now complete the sums so that the bounds are X]j=o • "^^^ crucial observations 
at this point are again that 4*-2*"^ = for i = -1, and likewise (^) 9*+ (^) 5' + 
(^) 3* = for i G {—1, —2}, and so the terms to be added are all 0. 

Using the binomial theorem to sum, we find 



J{A,k) = 16J(4,fc - 1) + (^yjlO'^ + [—j^'' + (6)4^^ + (8)^^- 

At this point the desired generating function can be shown to have coeflicients 
given by the same recursion. By iterating the recursion and summing the truncated 
geometric series, we obtain the formula given in Theorem 4. 

□ 

As expected, the values for the formula for J(4, k) are indeed zero at A: G 
{ — 1,-2,— 3}, and so a similar neat expansion holds for J(5,fc). The process of 
completing the binomial, summing with the binomial formula, and then iterating 
the recursion and summing the geometric series will provide explicit formulas easily 
as long as this property holds. 

From the expansion of Zn (t) at the end of the previous section it is clear that the 
sequences J{n, •) will all be tidy exponential sums with at most Fibonacci-many 
terms, possibly fewer as overlaps occur. If these sums are in fact zero in just exactly 
the necessary domain outside of the original meaningful domain of the indices, we 
think that this itself would be an interesting observation on the behavior of the 
combinatorics of the problem. 
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